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[57] ABSTRACT 

A technique based on synthetic aperture radar (SAR) 
interferometry is used to measure very small (1 cm or 
less) surface deformations with good resolution (10 m) 
over large areas (50 km). It can be used for accurate 
measurements of many geophysical phenomena, includ- 
ing swelling and buckling in fault zones, residual, verti- 
cal and lateral displacements from seismic events and 
prevolcanic swelling. Two SAR images are made of a 
scene by two spaced antennas and a difference interfer- 
ogram of the scene is made. After unwrapping phases of 
pixels of the difference interferogram, surface motion or 
deformation changes of the surface are observed. A 
second interferogram of the same scene is made from a 
different pair of images, at least one of which is made 
after some elapsed time. The second interferogram is 
then compared with the first interferogram to detect 
changes in line of sight position of pixels. By resolving 
line of sight observations into their vector components 
in other sets of interferograms along at least one other 
direction, lateral motions may be recovered in their 
entirety. Since in general, the SAR images are made 
from flight tracks that are separated, it is not possible to 
distinguish surface changes from the parallax caused by 
topography. However, a third image may be used to 
remove the topography and leave only the surface 
changes. 


13 Claims, 3 Drawing Sheets 
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METHOD FOR DETECT ING SURFACE MOTIONS 
AND MAPPING SMALL TERRESTRIAL OR 
PLANETARY SURFACE DEFORMATIONS WITH 

SYNTHETIC APERTURE RADAR 5 

ORIGIN OF THE INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA contract, and is sub- 
ject to the provisions of Public Law 96-517 (35 USC 
202) in which the Contractor has elected not to retain 
title. 

TECHNICAL FIELD 

15 

This invention relates to a method of detecting sur- 
face motions and measuring and mapping changes in 
terrestrial or planetary surface deformations with high 
resolution (1 cm) and high spatial resolution (typically 
10 m) over large areas (typically swaths 50 km wide) 20 
using synthetic aperture radar (SAR) images. The 
method is based on SAR interferometry and uses inter- 
ferograms produced from multiple images of the same 
area spaced apart in time. Changes in line of sight loca- 
tion are calculated from the measurement of the phase 25 
differences between corresponding pixels of the inter- 
ferograms. 

BACKGROUND ART 

In an airborne or spacebome SAR, microwaves illu- 30 
minate the terrain below the platform, and the echoes 
scattered from the surface are collected. Subsequent 
signal processing performed on the echoes produces an 
image of the scene, where each picture element (pixel) 
contains amplitude and phase information. A display of 35 
the amplitude information produces the familiar SAR 
images, which look very similar to photographs with 
high spatial resolution, typically — 10 m in a swath 50 
km wide. 

40 

It has been demonstrated that interferometric SAR 
observations provide a means for obtaining high resolu- 
tion topographic terrain maps from data acquired simul- 
taneously at two slightly displaced antennas. R. M. 
Goldstein, E. R. Caro and C. Wu, U.S. Pat. No. 4 _ 
4,551,724 issued Nov. 5, 1987; H. A. Zebker and R. M. 
Goldstein, “Topographic Mapping from Interferomet- 
ric Synthetic Aperture Radar Observations,” Journal of 
Geophysical Research, Vol. 91, No. B5, pp. 4993-4999, 
(1986). Typically, topographic maps are produced from 50 
a stereo pair of optical photographs in which the ap- 
pearance of the same terrain at slightly different posi- 
tions is used to deduce the height of the terrain. That is 
distinct from the interferometric approach used in the 
present invention. The following is excerpted from the 55 
paper cited above. 

In producing an interferometric SAR observation, 
two transversely spaced antennas Ai and A 2 are 
mounted on an aircraft equipped with side-looking 
radar as shown in FIG. la. Reflected energy from the 60 
ground is received simultaneously by the two antennas 
and separately recorded and processed to yield two 
~ 10-m resolution complex images (magnitude and 
phase). The phase <f> of the return signal from each 
resolution element (pixel) is maintained together with 65 
its magnitude. The two images are then combined pixel 
by pixel to obtain a single image with the phase of each 
pixel being the difference of the two original phases and 
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the magnitude of each pixel the product of the two 
original magnitudes. 

The radar system yields the line-of-sight distance 
(slant range) p and p' from each antenna to the target 
point P. The vertical distance h from the antenna to the 
surface is obtained from the solution of simple geomet- 
ric relationships using either slant range. However, the 
slant range is indeterminate within the slant range reso- 
lution (typically ±3.75 m) of the radar which depends 
principally on the transmitted signal bandwidth. The 
interferometric technique permits determination of the 
slant range p-p' to much greater precision (typically 1 
cm) as shown in FIG. lb. 

The two measurements of p and p', with the geomet- 
rical factors B and 0 shown in FIG. lb , are sufficient to 
find aircraft height h above each point in the radar 
image and from that deduce elevation of each point. 
The height h, when combined with the conventional 
along-track and slant-range (p-p') measurements, yields 
the three-dimensional position of a given point in the 
image. Knowledge of the location of each point then 
permits rectification of the image to a cartographically 
accurate form. In that manner, two SAR images re- 
corded simultaneously with two receiving antennas 
separated slightly in the cross-track direction, or re- 
corded with one antenna on two parallel but spaced 
apart passes, will provide topographic data of the ter- 
rain. 

As just noted, it is possible to perform the necessary 
measurements for SAR interferometry with only one 
antenna mounted on one orbiting spacecraft, or on one 
or two aircraft, by making two passes over the terrain 
along separate but parallel paths. If the two passes are 
not parallel and instead are skewed, the difficulty of 
detennining elevation of corresponding pixels in the 
separate images is substantially increased. A. K. Gabriel 
and R. M. Goldstein, “Crossed orbit interferometry: 
theory and experimental results from SIR-B,” Int. J. 
Remote Sensing, Vol. 9, No. 5, pp. 857-872 (1988). In 
any case, it is recognized that in determining the eleva- 
tion of any given point in the scene, the phase of the 
radar echoes may only be measured modulo 2tt, 
whereas the whole phase N277- is needed to determine 
absolute elevation. An approach to “unwrapping” the 
2tt ambiguities in the data set has been presented in 
detail. R. M. Goldstein, H. A. Zebker and C. L. Wer- 
ner, “Satellite radar interferometry: Two-dimensional 
phase unwrapping,” Radio Science, Vol. 23, No. 4, pp. 
713-720 (1988). The algorithm for deriving a satisfac- 
tory set of phases for SAR interferometry detailed there 
is incorporated here by reference. 

STATEMENT OF THE INVENTION 

A general object of this invention is to use multiple (2 
or more) SAR images for detecting surface motions, 
and mapping very small (~ 1 cm) surface deformations 
(changes in elevation) over very large areas. The ex- 
treme sensitivity of this technique to surface motion or 
deformation in terms of high spatial resolution (typi- 
cally 10 m) and broad swath coverage (typically 50 km 
wide) means it could be used to make extensive, accu- 
rate measurements of geophysical phenomena, includ- 
ing heaving and buckling in fault zones, tectonic plate 
motions, residual displacements from seismic events and 
motion from prevolcanic swelling. 

For detecting surface, two SAR images are mode of 
an area spaced in time, but with flight paths duplicated, 
such as by accurately duplicating the flight path later in 
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time or mounting a second SAR antenna in an aircraft 
well aft of the first antenna. After processing, a differ- 
ence interferogram of the area is made from two result- 
ing SAR images. Then, after unwrapping phases of 
pixels of the difference interferogram (to restore 2 tt 5 
multiples of phase that disappear in the making of the 
difference interferogram), the surface motion or defor- 
mation changes are observed. 

In the more general case, two interferograms are 
made from different pairs of three (or more) SAR im- 10 
ages of the same area (sometimes referred to hereinafter 
as the same “scene”)* where the second pair involves at 
least one image taken at a different time from images of 
the first pair. If the ground tracks for the SAR images 
are not parallel, an algorithm is used to compensate for 15 
this in the azimuth processing of data for the images, but 
before the interferograms are made, it is necessary to 
register the images. This is typically accomplished by 
resampling one image until it overlays the other accu- 
rately. Then an interferogram is made by calculating the 20 
phases corresponding to each pixel and computing their 
differences. To resolve phase ambiguity for precise 
elevation information, phase “unwrapping” is required. 

The unwrapped phases of the two interferograms are 
again differenced to produce a third interferogram re- 25 
ferred to hereinafter as a “double-difference interfero- 
gram.” This second differencing step removes phase 
changes from pixel to pixel due to topography, leaving 
a new phase image with nonzero phases only in areas 
where the terrain surface has moved in the time be- 30 
tween the observations for the first pair of images and 
the second pair of images used to produce the first and 
second interferograms. The phase differences of the 
double difference interferogram can be used to deduce 
the amount of surface deformation that has occurred. 35 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. la illustrates schematically the observational 
geometry of a SAR interferometer, and FIG. lb illus- 
trates the geometry when the baseline between antennas 40 
is not horizontal. 

FIG. 2a illustrates the geometry of two SAR obser- 
vations of a scene made from different locations in space 
using the same antenna for repeat observations from 
multiple passes of a satellite. 45 

FIG. 2b illustrates the geometry for measuring sur- 
face motion using two or three SAR observations of a 
scene where the surface motion is not along the line of 
sight. 

FIG. 3 illustrates the geometry of multiple observa- 50 
tion interferometry using three observations of a scene 
made from three locations separated by baselines p, q 
and r. 

FIG. 4 illustrates an ellipse tilted at 45° as determined 
from A<f>i2=(47r/X) p cos 0 and AcJ>i3=(4 tt/X) q cos 55 
(0-fa), where <f>=4 irp/\, and X is the radar system 
wavelength. Phase measurements used for double dif- 
ference interferometry are represented by the ellipse. 

FIG. 5 illustrates graphs used for interpretation of 
height changes of a site in the Imperial Valley of Cali- 60 
fomia. Watering is assumed to make the ground swell 
quickly over several hours followed by a slow shrink- 
age over several days. 

DETAILED DESCRIPTION OF THE 65 

INVENTION 

Referring to FIG. 2a, the range from two orbiting 
Seasat spacecraft antennas Ai and A 2 (spaced apart a 
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distance B perpendicular to the flight path) to a point P 
(scene resolution element) on the terrestrial surface is 
denoted p' for the first antenna Ai and p for the second 
antenna A 2 . The phases associated with these ranges are 

(f>—47Tp/X (1) 


where X is the radar system wavelength. 

For the present invention, these phases are taken to 
be the phases associated with the pixels in the two im- 
ages recorded before one interferogram is formed. This 
is not necessarily true since phase errors can occur in 
processing; however, as long as identical processing is 
performed on all points (pixels) of the images, the same 
errors are introduced everywhere, and constitute only a 
global phase constant that will disappear in the later 
step of phase differencing. 

In the following description of an example of the 
present invention, three Seasat L-band SAR images 
obtained of an area in the Imperial Valley of California 
over a period of 12 days were used to make two inter- 
ferograms, and then make a double-difference interfero- 
gram which clearly showed areas of surface deforma- 
tion changes. These changes were attributed to swelling 
of water-absorbing clays in the scene. 

In correlating the three images made from that raw 
data, it was necessary to process the data to be sure that 
azimuth focusing was identical for each of the three 
images, or errors would have been introduced into their 
phases with a resulting loss of phase coherence. This 
was accomplished by performing a conventional focus- 
ing algorithm on the first image, thereby to compensate 
for the radar looking forward or back (squinting); the 
resulting azimuth chirp rates were then used in the 
other two Seasat images as well. The processing for 
azimuth focusing should also include coregistration, 
radar motion compensation, and yaw, pitch, and roll of 
the platform, an application of the invention may re- 
quire. 

In order for the two sets of phases associated with a 
number of pixels to be mutually coherent, there are 
restrictions on the distance B between antennas. Specifi- 
cally, the distance B must be small enough that the 
speckle in the image from each pass does not change 
substantially; this requires both observations to be 
within the effective beamwidth of a reradiating resolu- 
tion element. This is discussed in the paper by Gabriel 
and Goldstein (1988) cited above. 

The Seasat L-band SAR had some orbits that met the 
criteria noted above for inteferometry. A full descrip- 
tion of the Seasat radar and correlator involved may be 
found in F. T. Ulaby, R. K. Moore and A. K. Fung, 
“Microwave Remote Sensing,” Addison Wesley, Read- 
ing, Mass., p. 48, 1981. 

If the orbits had overlapped exactly such that B=0, 
only two passes would have been necessary to detect 
surface motion or changes in surface deformation; how- 
ever, with a nonzero baseline, such changes become 
scrambled with the topographic information necessitat- 
ing at least a third SAR image. Sequential pairs of im- 
ages, such as images 1 and 2 for one pair and 2 and 3 for 
the second pair, are processed to form first and second 
single difference interferograms by multiplying the first 
image of a pair by the complex conjugate of the second 
image of the pair. 

This interferometric SAR technique measures only 
one-dimensional changes along the line of sight. If other 
sets of interferograms were made along at least one 
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other (nonparallel) direction, surface motions, including 
lateral motions, may be recovered in their entirety by 
resolving the observations into their vector components 
about predetermined reference axes that are mutually 
perpendicular with one of the three axes preferably 5 
selected to be vertical. This is illustrated in FIG. 2b 
which shows three different lines of sight LOSi, LOS 2 
and LOS 3 for interferograms used to measure one- 
dimensional changes in the position of a pixel from point 
Pi to point P 2 . By resolving the LOS components Ci, 10 
C 2 and C 3 of the motion in the direction of observation, 
total motion of the pixel to the new location P 2 from the 
original location Pi can be recovered from the geome- 
try knowing the angles involved which are measured 
and stored by the SAR system. Thus, two or three sets 15 
of interferograms may be used to measure surface mo- 
tion. The surface motion of a pixel is shown in FIG. 2b 
as being vertical for ease of illustration but any lateral 
motion may be recovered in the same manner once the 
observations of the pixel are resolved into vector com- 20 
ponents. 

Each of the phases from the original images is mea- 
sured to the nearest multiple of 27 r. Thus, two targets 
will appear at the same phase if their ranges differ by an 
integral number of wavelengths. Similarly, in an inter- 25 
ferogram, targets at different altitudes can appear at the 
same phase. It is necessary to find the 277 multiples of 
phase that disappeared or were “wrapped.” The proce- 
dure is referred to hereinafter as “unwrapping” and 
constitutes recovering the multiples of 27 t phase differ- 30 
ences that “disappear” in the interferogram. The gen- 
eral approach is to view the interferograms to check for 
fringes and then starting at some location (seed) in the 
interferogram where the phases have little noise, pro- 
ducing fringes that are bright and clear, integrating the 35 
adjacent pixel’s phase differences in expanding con- 
tours, and adding 2tt to the integrated pixel phases each 
time they become one more multiple of 2tt greater than 
the seed. Problems can arise where the image is noisy, 
or on steep slopes where the phase becomes noiselike, ^ 
because there are many resolution elements at the same 
range. This can lead to phase ambiguities for adjacent 
non-noisy areas; this problem is treated in the paper by 
Gabriel and Goldstein (1988) cited above. 

FIG. 3 illustrates the geometry associated with three 4g 
serial observations of a scene. The three observations 
are separated by baselines p, q and r. The phase associ- 
ated with a given resolution element (rezel) is given by 
Equation (1), and the phase difference obtained from 
the first two observations of that rezel is go 

A4>12= (4 tt/X)p cos d (2) 

Similarly, the first and third observations yield 

A4>i3-(4ir/\)^ cos (0+a) (3) 55 

Note that A ^ 12 /p and A(f>i 3 /q are the coordinates 
corresponding to the parametric equations of an ellipse 
FIG. 4. The positions of points PI and P2 on the ellipse 
are defined by sin a/(l— cosa) i and sin a/( 1 +cosa), 60 
respectively. In the case of the Seasat data used here, 0 
is confined to the narrow range 21 < 0 <23 and a = 180, 
so the ellipse almost degenerates into a straight line. If 
the orbits are not exactly parallel, the ellipse in FIG. 4 
will be slightly different at different locations in azi- 65 
muth. 

The two phases A^>n and A<f>i 3 cannot be differenced 
directly. First, each interferogram must be “un- 


6 

wrapped,” as described above. Then A<f>i 2 and AcJ>i 3 are 
plotted against each other, as in FIG. 4. Disturbances of 
the surface then show up as deviations of the points 
from the best-fit ellipse. 

Where the surface has been heavily disturbed 
(changes smaller than a pixel) as from plowing or a 
landslide, the phase sets have no relationship to each 
other. The resulting plot shows only randomly scat- 
tered points. However, where many rezels have been 
affected the same way (changes larger than a pixel), 
such as from a residual earthquake displacement, a regu- 
lar deviation from the ellipse will occur. The expression 

4> = A4m2/?-A<|>i3//>+c (4) 

is used to estimate the amount of phase deviation, where 
c is a constant, and p, q and c are taken from the best-fit 
ellipse for all of the data. Equation (4) highlights the 
aforementioned ambiguity which cannot be eliminated 
from only three observations. Motion of an area giving 
rise to changes in the phases of the associated pixels may 
have occured between observations 1 and 2 (thereby 
changing A<J>n) or may have occurred in the opposite 
direction between observations 1 and 3 (thereby chang- 
ing A<f>i 3 ). These two possibilities have different sensi- 
tivities that depend on the two baselines. One phase 
change is scaled by p/q, and the other by q/p. For the 
present example, this change was equal to 3 cm if it 
occurred during the long baseline observations, and 0.3 
cm if it occurred during the short baseline observations. 
Either a fourth observation and an obvious extension of 
the theory or some independent measurement will re- 
solve this ambiguity. 

The Seasat orbits in question were not exactly paral- 
lel, resulting in a very slight shearing of the images, 
which meant they could not be overlaid exactly. The 
second and third images were resampled to correct this 
using a procedure described in the paper by Gabriel and 
Goldstein (1988), cited above, and two single-difference 
interferograms were formed from the triplet. 

The phase of each interferogram was then un- 
wrapped and scaled by the baseline ratio as described in 
Equation (3), et seq., and the scaled interferogram pha- 
ses were differenced by multiplying the first single-dif- 
ference interferogram by the complex conjugate of said 
second single-difference interferogram, thereby pro- 
ducing a new image referred to as the ‘double-differ- 
ence interferogram.’ The slight divergence between 
orbits also produced a small phase gradient across this 
image and was corrected by fitting a plane to the result- 
ing phase difference map, which is the practical proce- 
dure for finding the small changes in the constants p, q, 
and c mentioned above. Thus, any linear phase errors in 
the along-track or the cross-track directions are re- 
moved. The consequence of this procedure is that only 
changes within a given image may be interpreted in 
terms of surface motion; in other words, comparing 
difference interferograms from two different sets of 
three orbits must allow for different phase gradients in 
each image. In the present example, the results were 
then interpreted according to known data about the 
imaging site. 

The Imperial Valley is an extremely busy agricultural 
area having many active fields that are constantly culti- 
vated. If a field was plowed or harvested, its surface 
changed and the speckle associated with the surface 
changed. Thus, two images of a field plowed between 
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observations yielded noise-like phase characteristics 
when they were interfered as the signals were no longer 
coherent. Such was the case with several of the fields in 
the experiment; their relative phases did not correlate. 
The fact that the noiselike phases follow field bound- 5 
aries is conclusive evidence that a physical change oc- 
curred; the effect is not from some systematic error in 
the experiment. Also apparent were fields with a shifted 
phase, which were the objects of interest, and large 
areas with no phase change, as would be expected. 10 

The site for the experiment was originally chosen 
because there were two small earthquakes (Richter 2.5) 
that occurred between the observations with epicenters 
near a radio tower. It may be safely concluded that any 
residual seismic displacements were less than the 1 cm 15 
accuracy of this technique. In the broader scene, there 
were clearly significant phase changes (90° phase 
shifts). Some physical change on the ground caused the 
phase to change in a coherent manner across entire 
fields. After some investigation, the most likely explana- 20 
tion was that those fields were watered between the 
observations. While the dielectric properties of the soil 
change with the addition of water, the Fresnel reflec- 
tion coefficient does not change enough to allow the 
large (90°) phase shifts seen unless there is penetration 25 
and subsurface scattering, in which case the large di- 
electric constant of water might cause a propagation 
delay through the wet soil. However, numerical esti- 
mates do not support this view. A better explanation is 
that the soils in the area are composed of hydroscopic 30 
clays that expand when wet. The phase changes ob- 
served probably resulted from expansion of an entire 
field when it was watered or contraction when it dried 
out. While such motions would not be large, the sensi- 
tivity of the method is so high that it could be visible. 35 

Each field served by the Imperial Valley Water Dis- 
trict has one or more gates that allow irrigation water 
from one of the surrounding canals. It was possible to 
collect some ground truth, consisting of ten year old 
watering records from the area. Records exist for ap- 40 
proximately 85% of the fields that exhibited a signifi- 
cant phase change. In addition, at each such site, re- 
cords were sought for at least two nearby fields that had 
no phase change. 

The watering records were interpreted using the 45 
model presented in FIG. 5, which shows the motion of 
the surface resulting from watering and indicates sche- 
matically the changes in propagation path associated 
with the motion. The water was assumed to be absorbed 
(causing swelling) over a period of a few hours and to 50 
evaporate (causing shrinkage) over a few days. Based 
on the watering records, which gave the date and 
amount of water supplied to each field, estimates of the 
amount of moisture present were made for each of 52 
sites in the image. Then Equation (4) was used to esti- 55 
mate the expected phase change, which was classified 
into seven categories from strongly negative through 
zero to strongly positive. Each of the 52 sites was simi- 
larly assigned to one of seven categories of phase 
changes, from negative to zero to positive. The two lists 60 
were then compared. In 48 cases, there was agreement 
between the theory and the experiment; in two cases, 
the result was ambiguous because the phase change was 
too noise for clear classification. In the final two cases, 
the phase change disagreed with what was predicted 65 
which was attributed to a clerical error in looking up 
the watering records. 

We claim: 
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1. A method for detecting surface motions or map- 
ping small terrestrial or planetary surface deformation 
changes over large areas using synthetic aperture radar 
comprising the steps of 

making two synthetic aperture radar images sequen- 
tial in time, but of the same area with flight paths 
substantially duplicated, 

processing said synthetic aperture radar images for 
azimuth focusing, co-registration, motion compen- 
sation, and yaw, pitch and roll, of the platform 
forming a difference interferogram of said area from 
said two synthetic aperture radar images, after 
processing in the previous step, 
unwrapping phases of pixels of said difference inter- 
ferogram by finding the 2 tt multiples of phase, if 
any, that disappear in the making of said difference 
interferogram, and 

observing line of sight surface motion or deformation 
changes in said terrestrial or planetary surface from 
said unwrapped difference interferogram, where 
said line of sight is of said synthetic aperture radar 
while making said images. 

2. A method as defined in claim 1 including the step 
of making at least one additional set of two or more 
sequential synthetic aperture radar images of the same 
area, but with flight paths producing different lines of 
sight for said synthetic aperture radar, recovering the 
line of sight component of said line of sight surface 
motion or deformation changes from each set of radar 
images, and from said line of sight components recover- 
ing total surface motion, including lateral surface mo- 
tion or surface deformation changes by resolving said 
line of sight components into vertical and lateral vector 
components. 

3. A method for mapping small terrestrial or plane- 
tary surface deformation changes over large areas with 
synthetic aperture radar comprising the steps of 

making three sequential synthetic aperture radar im- 
ages of the same scene with phase coherence be- 
tween corresponding pixels of said images, at least 
one of which is made significantly later than the 
other two 

processing said synthetic aperture radar images for 
azimuth focusing, co-registration, motion compen- 
sation, and yaw, pitch and roll, of the platform 
forming a first single-difference interferogram from a 
first pair of synthetic aperture images after process- 
ing in the second step, said first pair being selected 
from said three sequential synthetic radar images, 
forming a second single-difference interferogram 
from a second pair of synthetic aperture images 
after processing in the second step, said second pair 
being distinct from said first pair selected from said 
three sequential synthetic aperture radar images, 
unwrapping phases of pixels of said first and second 
single-difference interferograms by finding the 2 tt 
multiples of phase that disappear in the making of 
said first and second single-difference interfero- 
grams, 

forming a double-difference interferogram from said 
first single-difference interferogram and said sec- 
ond single-difference interferogram, thereby re- 
moving phase changes from corresponding pixels 
due to topography, leaving an image of nonzero 
phases only in areas where surface deformation has 
changed in the time between said first pair of im- 
ages and said second pair of images, and 



4,975,704 


9 

observing surface deformations in said terrestrial or 
planetary surface from said double-difference inter- 
ferogram. 

4. A method as defined in claim 3 wherein the step of 

forming said double-difference interferogram includes 5 
scaling fringes of said first and second single-difference 
interferograms by the ratio of baselines of the observa- 
tions used for making said three sequential synthetic 
aperture radar images for said first and second single- 
difference interferogram. 10 

5. A method as defined in claim 1 wherein said step of 
processing synthetic aperture radar images for azimuth 
focus includes resampling images until they overlay. 

6. A method as defined in claim 5 wherein said step of 
processing synthetic aperture radar images includes 15 
removing phase gradients caused by path divergence of 
said synthetic aperture radar during said observations 
used for making said sequential synthetic aperture radar 
images. 

7. A method as defined in claim 3 including the step 20 
of making at least one additional set of two or more 
sequential synthetic aperture radar images of the same 
area, but with flight paths producing different lines of 
sight for said synthetic aperture radar, recovering the 
line of sight component of said line of sight surface 25 
motion or deformation changes from each set of syn- 
thetic aperture radar images, and from said line of sight 
components recovering total surface motion, including 
lateral surface motion or surface deformation changes 
by resolving said line of sight components into vertical 30 
and lateral vector components. 

8. A method for mapping small terrestrial or plane- 
tary surface deformations over large areas using syn- 
thetic aperture radar comprising the steps of 

making multiple sequential synthetic aperture radar 35 
images of the same scene with phase coherence 
between corresponding pixels, 

processing said synthetic aperture radar images for 
azimuth focusing, 

forming a first single-difference interferogram from a 40 
first image of a first pair and a second image of said 
first pair, said first pair of images being selected 
from sequential ones of said multiple sequential 
synthetic radar images, 

forming a second single-difference interferogram 45 
from a first synthetic aperture radar image of a 
second pair and a second synthetic aperture radar 
image of said second pair, said second pair of syn- 
thetic aperture radar images being selected from 
sequential ones of said multiple sequential synthetic 50 
radar images, where each of said first and second 
synthetic aperture radar images of said second pair 
of synthetic aperture radar images is distinct from 
said first pair of synthetic aperture radar images, 
and at least one synthetic aperture radar image of 55 
said second pair is significantly later in sequence 


10 

than both synthetic aperture radar images of said 
first pair, 

unwrapping phases of each pixel of said first and 
second single-difference interferograms by finding 
the 2ir multiples of phase that disappear in the 
making of said first and second single-difference 
interferograms, 

forming a double-difference interferogram from said 
first single-difference interferogram and said sec- 
ond single-difference interferogram, thereby re- 
moving phase changes between corresponding 
pixels due to topography, leaving an image of non- 
zero phases only in areas where surface deforma- 
tion has changed in the time between said first pair 
of synthetic aperture radar images and said later 
observation, and 

observing surface deformation changes in said terres- 
trial or planetary surface from said double-differ- 
ence interferogram. 

9. A method as defined in claim 8 wherein the step of 
forming said double-difference interferogram includes 
scaling fringes of said first and second single-difference 
interferograms by the ratio of baselines of the observa- 
tions used for making said sequential synthetic aperture 
radar images for said first and second single-difference 
interferograms. 

10. A method as defined in claim 9 wherein said step 
of processing synthetic aperture radar images for azi- 
muth focus includes resampling synthetic aperture 
radar images until they overlay. 

11. A method as defined in claim 10 wherein said step 
of processing synthetic aperture radar images includes 
removing phase gradients caused by path divergence of 
said synthetic aperture radar during said observations 
used for making said sequential synthetic aperture radar 
images. 

12. A method as defined in claim 11 wherein the step 
of forming said multiple interferograms includes scaling 
fringes of said later interferograms by the ratio of base- 
lines of the observations used for making said sequential 
synthetic aperture radar images for said later interfero- 
grams. 

13. A method as defined in claim 8 including the step 
of making at least one additional set of two or more 
sequential synthetic aperture radar images of the same 
area, but with flight paths producing different lines of 
sight for said synthetic aperture radar, recovering the 
line of sight component of said line of sight surface 
motion or deformation changes from each set of syn- 
thetic aperture radar images, and from said line of sight 
components recovering total surface motion, including 
lateral surface motion or surface deformation changes 
by resolving said line of sight components into vertical 

and lateral vector components. 

* * * * * 
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